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SUMMARY 



Methods and charts are presented for computing sta- 
bility derivatives of a longitudinally straight V-bottom 
planing surface representing the forebody of a seaplane 
float or a flying-boat hull without chine flare. The 
charts for computing hydrodynamic derivatives were used 
in calculating the trim limit of stability for angles of 
dead rise of 10°, 20°, and 30°. The lower trim limit of 
stability of a seaplane planing on the forebody was cal- 
culated from measurements of lift, resistance, trimming 
moment, and wetted length of planing surfaces and. was 
found to be in good agreement with experimental values 
of the trim limit of stability* 

The velocity derivatives Z w and M w vera computed 
as a function of the draft, while the component due to 
the effect of vertical velocity on the trim was neglected. 
A comparison of the measured results with the calculated 
results indicated that the trim component was of minor 
importance and that better accuracy was obtained by neg- 
lecting it in the present calculations. 

INTRODUCTION 



The methods that are conventionally used in the anal- 
ysis of aerodynamic stability have been successfully 
employed in investigations of the stability of a seaplane 
in the planing condition. (See references 1 and 2.) The 
determination of both aerodynamic and hydrodynamic stabil- 
ity derivatives has been rather difficult.. In particular, 
the evaluation of the hydrodynamic components of the ve- 
locity derivatives (variations of lift and trimming moment 
with both linear and angular velocity) has involved assuniD 
tions of doubtful validity. These methods have given 



results that agree quAi^'Wv*if with the results of in- 
vestigations of the low-angle type of porpoising in which 
the afterbody is not involved. Applications of the 
methods to investigations of the high-angle type of por- 
poising appear more difficult and are not considered in 
the present report. 

The purpose of the present investigation was to 
compare values of the hydr odynami c derivatives as calcu- 
lated from general test data with the values measured 
from records of the disturbed motions of a planing "body 
and to develop a procedure for evaluating the derivatives 
more accurately than heretofore from general test data on 
planing surfaces. The four derivatives Z 2 , 2 V , Hg , 

and were measured from the damped oscillations of 

planing bodies. The derivative Z w was calculated and 
compared with the values obtained experiment ally % Of the 
four remaining hydr odynami c derivatives, Zg and M 2 
were obtained directly from general test data; Z^ and 

Viy were measured but not with sufficient accuracy for 
compari son with the calculati on a. as a check on the over- 
all accuracy of the methods and data, the lower trim 
limit of stability was computed for three different angles 
of dead ri se and was compared with the re suit s of t est s of 
dynamic models . 



The coefficient s used in the present report are 
defined as follows: 

/ \ 

Cl^ planing lift coefficient / 2 — \ 



SYMBOLS 





planing resistance coefficient 




C p center-of-pressure coefficient ( s/W . L 
C d draft coefficient (d/b) 



L.) 



C y speed coefficient (V/^/gb) 



C A load coefficient (A / wfc 3 ) 

C Aa initial load coefficient (A 0 /wb 3 ) 

where 

lp lift, pounds 

p density of water, slugs per cu"bic foot (1.97 slugs/ 
cu ft for water in NACA tank So. l) 

u horizontal velocity, feet per second 

"b team of hull, feet 

Ep resistance, pounds 

s distance from trailing edge to center of pressure, 
feet 

W.L. wetted length from trailing edge measured along keel, 
feet 

d draft , feet 

V speed, feet per second 

g acceleration of gravity, feet per second per second 

w specific weight of water, pounds per cubic foot 
(63.5 Ih/cu ft for water in 1TACA tank Ho. l) 

A load, pounds 

A 0 initial load, pounds 

Other syraools used in this report are defined as 
follows : 

k± horizontal distance of center of pressure forward of 
center of gravity 

k 2 vertical distance of center of pressure "below center 
of gravity 

I moment of inertia of seaplane alDout transverse axis 

through center of gravity, slug-feet 2 



mass of seaplane, slugs 

force along OZ-axi s , eq\ial to lift "but opposit e in 
sign 

trimming moment about transverse axis through center 
of gravity, positive when tending to raise "bow, 
pound-f cet 

period of oscillation of planing surface after 
disturbance, seconds 

angular velocity about center of gravity, radian s 
per second 

distance of perpendicular from center of gravity to 
keel forward of trailing edge, feet 

time, seconds 

resultant velocity of forebody due to small vertical 
velocity w impressed upon horizontal velocity u 

vortical velocity, feet per second 

vertical displacement from a condition of steady 
planing, positive downward , feet 

distance of pivot above water level 

, E n maximums on curve of disturbed motion of planing 
surface 

trim, radians 

angle of rotation about center of gravity, radians 



1 ogar it hmi c decrement 




1 dZ 

m d z 

i a_z 

m dw 



1 bZ 
m 3g 



M 



1 8M 
I dz 




1 dM 



I dw 




1 3M 

I ae 




1 dM 
I 3q 



CALCULATION OP DERIVATIVES 



Equations of Motion 



Perring and &lauert (reference l) demonstrated that 
the fore-and-aft degree of freedom may "be neglected in an 
analysis of porpoising and that the motion may "be treated 
as one having freedom in only trim and rise. In the con- 
ventional methods of stability analysis, the oscillations 
are assumed to he small. Maclaur in^s series is used to 
obtain expressions for the force and the moment and all 
derivatives of order higher than the first are neglected. 
The seaplane is assumed to be traveling at a constant 
forward speed under conditions of static equilibrium. An 
arbitrary disturbance of small amplitude is assumed and 
the equations of motion that describe the resulting small 
oscillations may be investigated either to determine the 
actual motions or to determine whether the motions diverge 
in amplitude or converge to zero. If the axes move with 
the seaplane and are taken with the origin at the initial 
position of the center of gravity, with OX forward and 
always parallel to the water surface and with OZ vertical 
and positive downward (fig. l), the equations of motion 
following a disturbance are 



m 



d % 



z 




36 



+ q 



3q 



(1) 



6 

I ^- = z — + w — + 6 — + q — ' (2) 

dt b z dw ae Sq. 

The division of equation (l) by m and of equation (2) 
"by I and the use of the notations 

Z z = 1 M 2 = i M ; and so forth 

m d z " I 9z 



give the equations for the acceleration 



3 
d s 

dt* 



d e 

dt 3 



= 2 Z Z + wZ w + 6Z g + qZ q (3) 



zll z + wM w + GMg + qM q (4) 



Basic Data for Computing Derivatives 

Data frcm reference 3 have "been plotted in a con- 
venient form that neglects the effect of Froudo^s number 
on lift, resistance, and the position of the center of 
pressure. Lift coefficient, resistance coefficient, and 
cent or-of~pressur e coefficient are plotted against draft 

coefficient; lift coefficient , resistan ce coefficient^ 

2 9 a 
draft coefficient draft coefficient 

and ce:iter-of~prossure coefficient are plotted against 
trim. Those plots are given in figures 2 to 19 as^faired 
curves for angles of dead rise of 10 , 20 , and 30 . 
Values for the draft wore computed from measurements of 
v/et tod lengths . 



Di spla cement Deri vat ives 

Deriva tive Z z The derivative Z s may he obtained 
from a plot of 0 L as a function of with T as 



7 



the parameter. The force Z is equal to the lift Lp 

tut opposite in sign and, at constant T and u, 
dz =5 dd. Hence, for unit mass, 



* 
i 

A 



7 - 1 91 



1 _ 
id d d 



1 ^ 



m BC d 



and Z z may he evaluated if the mass moving vertically 

is suhstituted in the equation and the slope ^ 
is ohtained from figure 2, 3, or 4. 

Derivative Zg ♦ - The derivative Zg is proportional 

to the rate of change of the force Z with change in 8 . 
A change in 8 implies a change in draft as well as in 
trim T ; hence , 

Z - A Si* - I (£-5. dd + 9Z dT\ 
6 m 66 m W d6 9t dG/ 

dr 

and, "because — L = 1, 
d6 

Zfi „ Z « . I^ilpuV 
b 2 d8 m dt 2 

In order to obtain — , which equals — , the equation 

d6 * dO 

relating d and 6 is differentiated. From figure 1, 

if 2 0 is the distance of the pivot above the water level, 

d a r sin 6 + p cos 6 - z 0 

where r is the distance of the pivot forward of the 
trailing edge and p is the distance of the pivot above 
the keel. Then, 



M = r cos 6 - P sin 6 
d& 



8 



and, for small angles, 

dd 

d8 = r ~ P9 

The sIotds - may be obtained from figure 5, 6, or 

7. Figures 5, 6, and 7 were obtained "by cress-plotting 
figures 2, 3, and 4, respectively. 

Trimming moments may "be plotted directly as a func- 
tion of draft or trim to obtain K z or Mg , but separate 
plots are required for each velocity to "be investigated* 
If the definitions of planing coefficient, resistance 
coefficient, and cent er-of-pre s sure coefficient are used 
ard if the effect of Froude l s law is neglected, these 
moment derivatives may he obtained at any speed from a 
plot of Op as a function of and t . 

Derivative M £ .- The expression for trimming moment 
i s — — — 

H = L p ki - R p K 2 

where k x is the horizontal distance of the center of 
pressure forward of the center of gravity and k 2 is the 
vertical distance of the center of pressure 'below the 
center of gravity. In order to obtain M z , this expres- 
sion is differentiated vlth respect to d, with T held 
constant : 

M i 3M m 1 A 3kx + ^ _ _ ka 

z I Bd I \ p Sd 1 Bd p dd 3d J 

IProm figure 1, 

ki = ( s - r) cos T + p sin T 

and 

k 2 s p cos T - (s - r) sin T 



9 



By dif f erentiation, 

akl c p ^ a a Qp 

dd re. T Bd 

and 

2£» s « _ d^£ 

dd 5 3d 

3R 

Tho partial derivative — £ may lie evaluated "by measur- 

3 d 

ing the slopes of the curves in figures 8 tc 10; likewise, 
— 2 is evaluated as a slope in figures 14 to 16. 

0 CL 

Derivative Mg The derivative Mg is analagous 
to Zq and may "be evaluated as the sum of two variables: 

6 ~ I 66 ~ I W 59 Bt d8' 

where 

9M Sk x BL p dkjg SR^ 

dr p 3t 6t p 3t dT 

Differentiating k x and k s with respect to T gives 



dk x dC p d ac p 
3t = T a * t 3t 



a * ~ TT - +* rt + p 



3ks BCp 
= - pt - dT — + r 

The -partial derivative may he evaluated from figures 

St 

11 to 13, and — E is obtained from figures 17 to 19. 

3t 



11 



dT 1 
dw u 



For unit mass and unit moment of inertia, the deriv- 
atives may te expressed in the form 



z w = LL z z + I z, 

U U 



u z u 

As will "be shown, in the following paragraphs, exper- 
imental data indicate that the components of Z w and M w 
which contain Zq and Mq are of minor importance and 
for the present purpose may "be neglected. 

Derivatives Zq and Kq.- In order to evaluate Zq 

and Mq, it is assumed that an angular velocity q may 
he considered to have the same effect upon the planing 
surface as the sum of the effects due to a vertical and a 
horizontal velocity - that is, 

1 / dZ du + dZ dw \ 
m \ 3u dq dw dq / 



= i az s 

* m dq 



Because 



I dq I ^ Bu dq dw dq ^ 



dq 



and 



12 



the two derivatives can be evaluated if ^— and ^- are 

3u 3u 

determined. 3y the use of 

L p = C Lp fpu 8 b 3 

— £ = C L _ pub 3 = — £ 

In order to obtain slopes are taken from a plot of 

du 

moment as a function of speed at constant load. In the 
expression 

M = k-Lp - k 0 R p 

L p is known, k ? . and k 3 havo already boen calculated, 
and ■ E p may be found by the use cf figures 8 to 10 because 
Cft is known. 

EVALUATION OP SOUTH' S DISCRIMINANT AND DETERMINATION 

01 CRITICAL TRIM 

If the hydrodynamic and aerodynamic derivatives have 
been found for a given load, speed, and trim, the deriva- 
tives are substituted in the equations cf motion. Routh's 
well-known criterions for the motion to be stable are 
that B, C, D, E, and H all be greater than zero, 
where 

B = -(Z w + M q ) 

C = -(M e + Z z - M q Z w + m w z 4 ) 
D = Z z M q - Z q M z + Z v Kg - Z 6 M W 
E = Z z M e - Z e M z 



and 



R = BCD - D -BE 



13 



The^lower limit of stability may be found when B , 
c » 3, and R have been evaluated for a number of 

angles . 



EXPERIMENTAL EVALUATION 0? DERIVATIVES 

The hydrodynamic derivatives 2 Z , Z w , Mq, and Mq 

were measured "by the "froe-oscillation method 11 that is 
sometimes used in measuring aerodynamic derivatives. A 
model composed of a planing surface and a tail plane (fig. 
21), which represent the forebcdy and the horizontal tail 
of a seaplane (described in reference 4), was first towed 
free to move vertically and locked in trim and vas then 
towed free to trim but locked in draft. In each case, 
the model was disturbed momentarily and a time history of 
the subsequent oscillations was obtained. If the oscilla- 
tions are assumed to be those of a damped linear oscilla- 
tor, the equation of motion with the model locked in trim 
is 

m = Z w — + Z 3 2 

dt s dt 

wh e r e 

Z w damping factor (2m 



Z 2 displacement factor 
and 

3? period of one oscillation 

6 logarithmic decrement (log 

z lf z 2 two consecutive maximums on curve 

A typical trace of the damped motion in draft with locked 
trim is given in figure 22. Analagous equations apply 
when the model is free to trim but locked in draft. 



• 6 -) 



m 



dt dt 



14 



where 



and 



M . 215 
P 



m 6 = i 



(4tt 2 + S 3 ) 

P ai ' 



A limited number of records were taken to obtain these 
four derivatives and the results are compared with the 
corresponding calculated values of the derivatives in 
tables I and II. The results given in table II show good 
agreement between the measured values of Z w and the 

value calculated by the formula 

m u 

The agreement is better and the calculations are simpler 

than if the term' — Z c is added. 

a u c 

LOW38E THIM LIMIT 03P STABILITY 



The methods and charts for computing the derivatives 
were applied to the specific problem of calculating the 
lower "trim limit of stability for angles of dead rise of 
10°, 20°, and 30°. The calculations were generally sim- 
ilar to those of reference 4 but were more extensive and 
employed the method of calculating the velocity deriva- 
tives previously described herein. The calculations were 
made for a model with a horizontal tail plane and without 
a wing. The dimensions, the mass, the moment of inertia, 
and other characteristics assumed for the model were the 
same as those of reference 4 and were as follows: 



15 



Mass, slugs 5.05 

Moment of inertia, slug-feet 3 5.2 

Beam, feet 1 .33 

Center of gravity, fraction bean above keel .... 1.25 

Center of gravity, fraction beam forward of T.3.. . .33 

Tail area, square feet. 3.47 

Aspect ratio of. tail plane 3.4 

Tail arji (measured from eg* to quart er- 

chord point), feet 3.93 

El e vat or area , square inches 492 

Elevator chord, percent total chord 48 



Table III presents the results of the calculations 
for one set of conditions: namely, angle of dead rise of 
30°, speed of 40 feet per second, and load of 60 pounds. 
The critical trim , - that is, the trim at which the calcu- 
lations indicate a transition from stability to instabil- 
ity, was determined at two speeds for four loads and for 
two angles of dead rise, 20° and 30°. The calculated 
values are plotted in figure 2Z. Graphs showing all the 
hydr o dynamic derivatives used in these comput ations are 
available, on request, from the Nat i on-al Advi sory C ommit t ee 
for Aeronautics. The faired curves of figure 2-3 are from 
reference 4 and were determined experimentally. The cal- 
culated values agree with the experimental values within 
the probable limits of experimental error, except at a 
speed of 30 feet per second and a load of 1 00 pounds . The 
di screpanoy for t hi s one point may be due to the fact that 
the positions of the center of pressure were obtained by 
extrapolating the data in reference 3 from 80 pounds to 
100 pounds. Also, the assumption that Proude's law of 
comparison may be neglected probably introduces significant 
errors at the low speeds and heavy loads. 

CONCLUDING REMARKS 



Methods and charts are presented herein for comput- 
ing the stability derivatives of a longitudinally straight 
V~bottom planing surface representing the forebody of a 
seaplane float or a flying-boat hull without chine flare. 
The methods are believed to be generally applicable in 
computing stability derivatives of a single planing sur- 
face and the charts for angles of dead rise of 10°, 20°, 
and 30° are believed to be satisfactory approximations 
for calculations of the lower trim limit of stability 
when the straight portion of the forebody forward of the 
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step is the only planing area involved. In particular, 
the results indicate that the velocity derivatives may he 
calculated satisfactorily as functions of the draft with 
trim constant while the component due to the effect of 
vertical velocity on the trim is neglected. 



Langley Memorial Aeronautical Lahoratory, 

National Advisory Committee for Aeronautics, 
Langley Pield, Va. 
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Figure 7.- Variation of j ^ a ^ n S coefficie nt witb t ^ Qne curye ?Qr 
e draft coefficient^ 

all values of Cd« Angle of dead rise, 30°. 
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draft coefficients 
Angle of dead rise, 10°. 
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Figure 12.- Variation of ™ ? 1 s tanc iL.coef £i c I ant with trim< 

draft coefficient2 
Angle of dead rise, 20°. 
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Figure 13.- Variation of resistance coefficient with trim< 

draft coefficient2 
Angle of dead rise, 30°. 
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Figure 15.- Variation of center-of -pressure coefficient with 

draft coefficient. Angle of dead rise, 20°. 
C p = s/W.L. > Cd - d/t>. 
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Figure 16.- Variation of center- of-pres sure coefficient with draft 
coefficient. Angle of dead rise, 30°. 0. Q = s/W.L. ; 

c d . = d/t). y 
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Variation of center-of -pressure coefficient with 
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Figure 20.- Resultant draft and trim of a planing forebody due 

to a small vertical velocity w impressed upon the 
horizontal velocity u. 
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Figure 23.- Comparison of calculated critical trim with experi- 
mental critical trim. Model with tail. 



